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ABSTRACT 

An  acoustic  energy  reflection  experiment  was  conducted  in  a  labora- 
tory tank  containing  water  and  an  artificial  kaolin  sediment.   The 
reflection  coefficient  for  the  sediment  was  measured  for  a  limited 
range  of  angles  of  incidence  at  one  frequency  and  one  set  of  measured 
sediment  mass  physical  and  viscoelastic  properties.   A  mathematical 
reflection  model  was  derived  in  terms  of  the  known  sediment  properties. 
The  measured  and  the  theoretically  calculated  reflection  coefficients 
are  compared.   The  measured  reflection  coefficient  and  a  probable  angle 
of  intromission  showed  poor  agreement  with  theoretically  predicted 
values.   For  the  sediment  parameters  of  the  relatively  soft  clay  used 
here,  the  reflection  model  is  quite  sensitive  to  values  of  sediment 
density  and  sound  speed  and  is  only  mildly  sensitive  to  values  of 
dynamic  rigidity  and  sound  absorption  coefficient  near  the  angle  of 
intromission. 
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I.   INTRODUCTION 

A.   BACKGROUND 

The  development  and  use  of  bottom  bounce  sonar  and  long  range 
passive  detection  systems,  and  the  use  of  sonar  in  shallow  water  all 
require  knowledge  of  the  effect  the  ocean  bottom  has  on  sound  propa- 
gation loss.   This  effect  is  most  often  characterized  by  the  quantity 
bottom  loss  (BL)  in  decibels,  which  is  related  to  the  reflection  co- 
efficient (R)  ,  the  ratio  of  the  reflected  acoustic  intensity  to  the 
incident  acoustic  intensity,  by  the  equation  BL  =  -10  log    R.   The 
large  variations  in  bottom  loss  throughout  the  world  ocean  bottoms, 
coupled  with  the  passive  or  active  sonar  figure  of  merit  sensitivity 
to  small  changes  in  this  parameter,  dictate  a  need  for  accurate 
measurement  of  bottom  loss  characteristics  in  each  area  of  the  oceans 
where  bottom  sediment  parameters  are  different. 

A  direct  method  for  collecting  such  information  would  be  an 
oceanographic  survey,  of  massive  proportions,  to  determine  the  acoustic 
reflection  coefficient  at  various  angles  of  incidence  and  various 
frequencies  of  interest  in  each  distinct  bottom  region  of  the  world's 
oceans.   On  the  other  hand,  surveys  of  other  physical  characteristics 
of  the  ocean  bottoms  have  been  going  on  for  years,  with  the  result 
that  some  core  samples  and  raw  and  analyzed  data  exist  for  practically 
every  oceanic  bottom  region.   Although  much  information  is  available 
in  the  literature  concerning  the  relation  of  reflection  and  transmission 
of  acoustic  energy  to  the  properties  of  the  reflecting  and  transmitting 


media,  a  round  table  of  exploration  geophysicists  concluded  recently 
that  their  understanding  of  the  "nature  of  seismic  reflection  is  far 
from  complete."*   In  anti-submarine  and  submarine  warfare,  where  the 
exact  nature  of  the  reflection  and  transmission  of  acoustic  energy 
at  the  ocean  bottom  must  be  well  known,  models  for  calculation  of  re- 
flection coefficients  generally  are  limited  to  consideration  of  a  plane, 
homogeneous,  compressional  wave,  transmitted  at  various  angles  of  in- 
cidence in  water  and  reflected  from  and  refracted  in  layered  sediments 
having  various  densities  and  other  physical  properties,  including 
acoustic  energy  attenuation. 

The  energy  loss  a  plane  compressional  wave  undergoes  upon  reflection 
from  the  ocean's  bottom  is  examined  in  a  number  of  reports,  papers  and 
textbooks.   Three  of  these  references  were  chosen  to  describe  the  nature 
of  the  problem  examined  in  this  paper.   White  [1]  was  used  as  a  basic 
text  which  describes  the  energy  loss  as  a  result  of  both  compressional 
and  shear  waves  being  refracted  in  the  sediment.   Hamilton  [2]  uses  the 
physical  properties  of  density,  compressional  and  shear  wave  speeds, 
along  with  related  values  of  bulk  modulus  (K)  and  modulus  of  dynamic 
rigidity  (y)  in  the  sediment  to  predict  values  of  reflection  coefficient 
and  bottom  loss.   Hamilton  treats  the  sediment  as  both  an  elastic  and 
a  viscoelastic  medium.  Bucker  et  al [3] compares  measured  and  theoretical- 
ly calculated  values  of  bottom  loss  for  multi-layered  viscoelastic 
ocean  sediments. 


* 
Anstey,  N.  A.,  and  Allen,  S.  J.,  "What  Direction  Should  We  Set 
for  Hydrocarbon  Exploration?"   Exploration  Geophysics,  Vol.  21,  No.  3, 
p.  412-421,  September  1973 


A  number  of  the  acoustic  reflection  experiments  that  have  been 
conducted  at  sea,  such  as  those  referenced  in  [2]  and  [3] ,  have  shown 
varying  degrees  of  agreement  between  theory  and  experiment.   The  dis- 
agreement usually  is  greatest  near  the  critical  angle  or  angle  of  in- 
tromission.  The  theoretical  modeling  of  these  experiments  has  been 
complicated  by  sediment  layering  and  other  non-uniformities  and  by 
difficulty  in  obtaining  realistic  physical  properties.   A  laboratory 
experiment  involving  a  smooth-surfaced  homogeneous  sediment  and  an 
undisturbed  water  layer  helps  to  eliminate  some  of  these  complications 
in  comparing  measured  and  theoretical  reflection  coefficients. 

B.   DYNAMIC  RIGIDITY 

A  key  element  in  improved  acoustic  reflection  models  has  been  the 
inclusion  of  the  dynamic  rigidity  and  dissipation  properties  of  the 
sediment.   Most  researchers,  such  as  Hamilton  [2]  and  Bucker  et  al  [3] , 
used  measured  Stoneley  wave  speeds  to  calculate  the  real  part  of  dynamic 
rigidity.   Only  estimates  of  the  imaginary  part  were  available.   Also, 
Stoke 's  assumption  that  in  a  viscoelastic  fluid  the  bulk  viscosity  is 
equal  to  zero  was  employed  by  these  researchers  to  estimate  a  value  of 
the  imaginary  part  of  the  first  Lame  constant  (A) . 

At  the  Naval  Postgraduate  School,  the  torsional  oscillator  method 
has  been  used  to  measure  the  real  and  imaginary  parts  of  the  complex 
dynamic  rigidity.  Cohen  [4]  described  one  set  of  these  measurements 
for  a  saturated,  water-kaolin  artificial  sediment,  and  Martinek  [5] , 
conducting  research  on  the  same  type  of  sediment,  measured  values  of 
absorption  coefficient. 


C.   REFLECTION  COEFFICIENT 

An  objective  of  this  research  was  to  examine  the  usefulness  and 
validity  of  the  values  of  dynamic  rigidity  and  absorption  coefficient, 
along  with  other  data  on  the  physical  properties  of  the  sediment,  in 
calculations  with  a  viscoelastic  reflection  model,  and  to  compare  the 
calculated  results  with  experimental  reflection  coefficients.   The 
viscoelastic  reflection  model  of  Bucker  et  al  [3]  was  modified  to  avoid 
the  use  of  the  assumption  that  the  bulk  viscosity  is  zero  for  a  compres- 
sional  wave.   Some  preliminary  experimental  measurements  were  made  of 
the  reflection  coefficient  at  one  frequency  and  for  one  sediment  condi- 
tion for  comparison  purposes. 
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II.   THEORETICAL  REFLECTION  COEFFICIENT 

A.   DESCRIPTION  OF  THE  MATHEMATICAL  MODEL 

Figure  1  illustrates  a  plane  wave  reflection  model  which  was 
assumed  to  be  an  approximation  to  the  experiment.   The  symbols  and 
notation  defined  in  Figure  1  were  drawn  from  Bucker  et  al  [3] .   Refer- 
ring to  this  figure,  (J)  and  ^  are  potential  functions  defined  in  the 
following  paragraph  and  x  is  positive  to  the  right  while  z  is  positive 
downward . 

A  summary  of  the  development  of  water-sediment  interface  equations 
in  [3]  follows.   The  wave  equation  is  given  by  the  relationship 
V  4>  +  (— )   (J)  =  0,  and  its  assumed  potential  function  solutions  are:  (1) 

4>'2  =  A'   exp  i(kx  +  a^z  -  cot),  (2) 

(J)1   =  A'   exp  i  (kx  +  a  kz  -  cot)  ,  (3) 

<J>"2  =  A"2  exp  i(kx  -  a2kz  -  wt)  ,  (4) 

V   =  B'   exp  i(kx  +  b  kz  -  cot)  (5) 

II  l        / 

A"1  -  0  ,     B"   =  0  . 

In  water  only  the  real  values  of  k  are  considered  for 

co  sin  G  . 

k  = .  (6) 

c 
w 

The  term  c   is  the  measured  speed  of  a  compressional  wave  in  water, 
A'   and  A"   are  the  measured  amplitudes  of  the  incident  and  reflected 
waves,  respectively,  A1   and  B'   are  unknown  quantities,  6.  is  the  angle 
of  incidence,  and  <f>"   and  H7"  are  both  equal  zero. 
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Figure  1.   Theoretical  Reflection  Model. 
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For  the  water-sediment  interface  two  boundary  conditions  may  be 
applied  to  the  potential  function: 

1.  the  vertical  particle  velocities  (w)  are  continuous  across  the 
interface; 

W  -  -^   +    -r-  (7) 

dz  dx 

2.  the  vertical  stress-tensor  elements  are  continuous  across  the 
interface; 

P    =  i(-P0)2(J)  -2U  ^  +  2)1  ^fl_)  (8) 

ZZ    W  dx2       8x82 

p   =  0  . 
zx 

Following  Bucker  [3] ,  three  interface  equations  are  derived  as  follows 
from  the  boundary  conditions  and  treatment  of  the  sediment  as  a  visco- 
elastic  half-space; 


a2(A'2  -  A"2)  -  a1  A^  +  B^  (9) 

P2  (f)2  (A'2  +  A"2)  =  yi[(b12-l)  (A^)  +2b1(B'1)]        (10) 

2ai  (A^)  =  (b^-1)  B^  (11) 


where: 


2 

2     Pl  C 

•i  ■    v^"1 

2 

2      Pl  ° 
b/  =    — -  1  (13) 

1        ^1 
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c 

c  =  _ —  -  horizontal  phase  speed  in  water 

sm  8 . 
1 


a  =  cot  9 
2        i 


p   =  density  of  sediment 
p  =  density  of  water 


The  Lame  constants  A  and  y  are  complex, 


Al  =  A'l  "  iA"l  U4) 

Pl  =  y'l  '  iy"l  "  U5) 


B.   SEDIMENT  PARAMETERS 

The  above  model  permits  a  solution  to  the  wave  equation  in  terms  of 
sediment  parameters  A,  y,  p  and  c,  of  which  A  is  the  only  term  which 
does  not  have  independently  determined  values.   Cohen  [4]  measured 
values  of  dynamic  rigidity,  Ji'   and  y"   for  a  range  of  sediment  densi- 
ties.  Martinek  [5]  measured  the  absorption  coefficient,  a  ,  and  the 
compressional  wave  speed  in  the  sediment,  c  ,  at  a  sediment  density 
comparable  to  Cohen's  and  at  frequencies  in  the  range  of  experimental 
interest.   These  values  are  summarized  in  Table  I  of  Section  IV.   Addi- 
tional sound  speed  and  density  values  were  measured  as  part  of  this 
experiment  and  recorded  in  Table  II  of  Section  IV. 

In  the  absence  of  scattering  from  the  water-sediment  interface,  a 

P 

in  theory  is  related  to  A"   and  y"   and  c   is  related  to  A '   and  y'  . 

With  these  known  quantities,  the  next  series  of  steps  transposes  the 

equations  to  eliminate  A'   and  A"   as  independent  variables  and  solves 

for  a  reflection  coefficient  in  terms  of  sediment  parameters  y1  ,  y"  , 

p  ,  c  ,  along  with  c  ,  po  and  0.  in  the  water  layer. 
-LI  w    2      i  J 
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C.   CALCULATION  OF  REFLECTION  COEFFICIENT 

In  the  following  development,  the  symbols  are  similar  to  those  used 
by  Borcherdt  [6]  but  with  a  reversed  sign  convention  for  the  complex 
exponential.   The  sign  convention  used  here  is  the  same  used  by  Bucker 
et  al  [3],   Definition  of  terms  is  as  follows: 

k  =  k1  -  ik"  =  —  =  complex  propagation  constant 
a  =  a'  -  ia"  =  complex  wave  speed  and 
CO  =  angular  frequency 


a  =  p.v, 


Kl  +  3  ^1 


1/2 


where  K  is  the  bulk  modulus  and  p.v.  indicates  principal  value. 

4 
From  Hamilton  [2]  ,  K  +  —  y,  =  A   +  2u    so  that 


a  =  p„v. 


A1  +  2y1  \       1/2 


(16) 


The  compressional  wave  speed  in  the  sediment,  c  ,  is  given  by 


-  =  Re   [  £■  ]   . 

c         a 


(17) 


is 


The  compressional  wave  absorption  coefficient  in  the  sediment 


a  =  calm  [  —  ] 
p        a 


(18) 


The   complex  conrpressional   wave   propagation   constant    is 


coa'  .   ooa"         co 

k   =  -] — T2   +    i  -i — T2    =  —  +    ia 

hi  H         c±  P 
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For  physical  reasons,  a   >  0  to  ensure  decreasing  values  of  the 
potential  function  with  positive  increasing  values  of  distance  along 
the  propagation  path.   By  the  same  physical  reasoning,  c   >  0. 

Algebraic  steps  to  solve  for  reflection  coefficient,  R,  in  terms  of 
u '  ,  U",  c  ,  a   and  frequency,  f,  and  independent  of  A*   and  A"   are 
shown  in  detail  in  Ax^pendix  A.   These  steps  yield  the  following  set 
of  equations: 


Z  = 


R 


A'   +  A" 
2 2 

A1   -  A" 
A  2      2 


A»2 

2 

Z-l 
Z+l 

2a2yik 


2     2 
P2(bl  +1)W 


(b2-D2 
^ +  2bi 


\ 


(A-10) (A-ll) 


H  +  iJ 


(A-9) 


(A-12) 


F  +  iG 


(A-15) 


F  = 


2  , 
P1C  U  1 


(y»1)2+(y"1)2 


-  l 


(A-16) 


Plc2y"l 


(y'1)2+(y"1)2 


(A-17) 


Pxc 


2/  1 


-  i<^>2 


(A-18) 


c  a 

(JLe, 


c  a   2 
l-(-iP) 


H  +  1 


(A-19) 
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Equation  (A-ll)  for  R  then  was  solved  using  standard  WATFORG  on 
the  IBM  360  computer  at  the  Naval  Postgraduate  School.   Values  of  6. 
were  chosen  from  zero  to  tt/2  and  the  physical  properties  of  the 
sediment  were  taken  from  Cohen  [4]  and  Martinek  [5] .   Values  for  BL 
and  R  were  computed  over  the  range  of  experimental  error  measured  in 
the  physical  properties  cited  in  these  two  references,  in  order  to 
examine  them  for  sensitivity  to  changes  in  these  values. 
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III.   MEASUREMENT  OF  REFLECTION  COEFFICIENT 

A.  SUMMARY  OF  THE  EXPERIMENT 

The  experiment  was  conducted  by  reflecting  a  pulsed  beam  of  sound 
from  a  water-sediment  interface  at  various  angles  of  incidence.   The 
incident  and  reflected  wave  amplitudes  were  detected  with  a  directional 
receiver  of  design  identical  to  that  of  the  transmitter.   A  sound  fre- 
quency of  120  kHz  was  used.   The  received  signals  were  amplified,  fil- 
tered and  signal  processed  to  improve  signal-to-noise  ratio.   The 
reflection  coefficient  was  computed  from  the  ratio  of  the  squares  of 
the  voltage  amplitude  of  the  direct  and  sediment  reflected  pulses, 
using  appropriate  corrections  for  spherical  divergence. 

B.  THE  APPARATUS 

The  artificial  sediment,  a  commercially  available  pure  kaolinite, 
was  contained  in  a  large,  rectangular  tank,  2.42  m  long  by  1.20  m  wide 
by  1.00  m  deep.   The  leveled,  saturated  sediment  thickness  was  60  cm. 
The  water  depth  was  31.5  cm. 

Identical  transmitting  and  receiving  transducers  were  set  up  in 
the  water  and  mounted  on  a  horizontal  rail  above  the  tank.   Each  trans- 
ducer could  be  rotated  through  360°  independently  in  both  the  vertical 
and  horizontal  axes  and  moved  horizontally  and  vertically  in  the  tank. 
The  two  electrostatic  type  transducers,  built  by  research  physicist 
Donald  F.  Spiel  of  the  Naval  Postgraduate  School,  have  an  active  circu- 
lar face  of  6.0  cm  diameter,  consisting  of  an  aluminum  coated,  one  mil 
thick  plastic  film  mounted  over  an  aluminum  electrode  which  is  set  in 
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a  larger  plastic  disc.   In  order  to  keep  the  air  gap  between  the  plastic 
film  and  the  electrode  uniform,  equalization  of  pressure  between  this  gap 
and  an  air  filled  space  behind  the  electrode  was  achieved  by  small  holes 
drilled  in  the  aluminum  electrode.   The  inside  of  the  transducer  was 
partially  evacuated.   At  120  kHz  the  3  dB  beam  width  is  about  20°. 

The  basic  geometry  of  the  experiment  is  shown  in  Figure  2.   Referring 
to  this  figure,  both  the  transmitting  and  receiving  transducers  were 
rotated  to  a  selected  angle  of  incidence,  0.. 

Based  on  the  choice  of  centerline  height  of  each  transducer  being 
23.5  cm  above  the  sediment  layer,  the  idealized  ray  path  y   was  computed 
for  each  0   and  the  transducers  were  positioned  along  the  horizontal 
rail  so  that  y  =  y  and  x   -  x  .   The  horizontal  placement  also  was 
chosen  to  ensure  that  the  same  region  of  sediment  was  ensonified  for  each 
case. 

The  range  of  angles  accessible  for  this  experiment  is  limited  by 
the  dimensions  of  the  apparatus  and  the  necessity  to  keep  the  sediment 
reflecting  surface  in  the  far  field  region  of  the  directional  trans- 
ducers.  The  extent  of  the  far  field  region  was  determined  experi- 
mentally by  plotting  the  directly  propagated  sound  field  pressure  as 
a  function  of  range. 

The  directivity  of  both  transducers  at  120  kHz  provided  a  suffi- 
cient reduction  in  off-axis  signals  to  eliminate  interference  from 
side  wall  and  water-air  surface  reflections.   This  was  confirmed  experi- 
mentally.  At  angles  above  70°  a  small  amount  of  interference  between 
direct  and  sediment  reflected  pulses  could  be  observed,  the  effects  of 
which  are  discussed  later." 
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Figure  2.   Experiment  Geometry. 
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C.  INSTRUMENTATION 

Figure  3  illustrates  the  transmitting  and  receiving  equipment  and 
circuitry  in  block  diagram  form.   This  configuration  is  similar  to 
that  used  by  Martinek  [5] .   The  Hewlett  Packard  204C  oscillator  fre- 
quency was  set  at  120  kHz  within  +  3  Hz  and  was  monitored  with  a  fre- 
quency counter.   The  General  Radio  1396-A  tone  burst  generator  was 
set  to  produce  two-cycle  tone  burst  pulses  with  an  approximate  10  msec 
delay  between  pulses.   This  delay  was  sufficient  to  prevent  undesired 
reflections  from  disturbing  the  measurements.   The  Hewlett  Packard 
467A  amplifier  was  used  to  boost  the  voltage  by  a  factor  of  ten  before 
the  signal  went  to  a  1:3  transformer  for  further  amplification.   The 
polarizing  circuits  for  both  the  electrostatic  transducers  are  shown 
in  Figure  4. 

The  received  signal  was  amplified  either  20  or  40  dB  by  a  Hewlett 
Packard  465A  amplifier  and  then  passed  through  Krohn-Hite  high  and  low 
pass  filters  set  to  a  20  kHz  pass  band  centered  at  120  kHz.   A  Prince- 
ton Applied  Research  (PAR)  model  160  boxcar  integrator  was  used  for 
signal  processing.   The  time  delay  feature  of  the  boxcar  integrator 
was  used  to  determine  the  times  of  pulse  arrival.   The  enhanced  signal 
from  the  arriving  pulse  was  then  fed  to  a  Moseley  X-Y  Recorder.   Fig- 
ure 5  which  was  traced  from  original  records  represents  typical  received 
wave  forms  after  integration  and  recording 0   Wave  form  points  selected 
to  determine  zero  to  peak  voltages  are  noted  on  this  figure. 

D.  DATA  COLLECTION  PROCESS  AND  ERROR  ANALYSIS 

A  number  of  angles  of  incidence  were  selected  and  appropriate  hori- 
zontal separations  were  pre-calculated.   At  each  setting  of  horizontal 
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spacing,  measurements  were  recorded  of  the  received  pulses  for  both  a 

direct  path  between  source  and  receiver  (6.  =0)  and  for  the  sediment 

1 

reflected  path.   Additional  gain  was  required  for  the  latter  measure- 
ment.  The  receiver  then  was  moved  its  full  horizontal  travel  to 
determine  that  no  other  pulses  were  being  reflected  from  below  the 
water-sediment  interface. 

Referring  to  Figure  5,  zero  to  peak  voltage  was  determined  by 
averaging  the  maximum  positive  and  negative;  peaks  of  each  recorded 
pulse.  In  all  cases  these  peaks  occurred  in  the  third  cycle  of  the 
pulse.   Voltage  output  anomalies  in  the  tone  burst  generator,  trans- 
formers, filters,  amplifiers,  boxcar  integrator  and  the  X-Y  recorder 
were  assumed  to  be  the  same  for  both  direct  path  and  reflected  path 
voltages  and  therefore  were  cancelled  by  division  in  the  reflection 
coefficient  ratio.   Also,  because  of  this  ratio  computation,  trans- 
ducer sensitivities,  sound  pressure  level  calculations,  etc.,  were 
not  required. 

The  log-log  plot  of  direct  path  average  zero  to  peak  amplitude  ver- 
sus horizontal  separation  was  constructed  to  determine  the  minimum 
horizontal  separation  at  which  the  pressure  varies  linearly  with  the 
inverse  range  (1/r) .   This  established  the  inner  limits  of  horizontal 
separation  which  in  turn  set  the  minimum  angle  of  incidence  at 
0.  =  50°.   Since  the  reflected  path  was  longer  than  the  direct  path, 
this  separation  insured  that  the  ensonified  patch  of  sediment  was  also 
in  the  far  field.   The  dimensions  of  the  tank  limited  the  maximum 
angle  of  incidence  to  6 .  =  76°. 

The  arrival  time  of  the  sediment  reflected  pulse  occurred  later 
than  the  arrival  time  of  the  direct  pulse.   The  measured  arrival  time 
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differences  between  the  two  pulses  were  compared  with  time  differences 
calculated  using  the  idealized  geometry  ray  path  difference  and  the 
known  speed  of  sound.   Compatibility  between  observed  and  calculated 
delays  confirmed  detection  of  the  sediment  reflected  signal. 

The  uncertainty  in  the  angle  of  incidence  (6.)  on  each  transducer 
is  about  +  1.0°.   Part  of  this  is  due  to  the  precision  with  which  the 
scale  could  be  read  and  part  arises  in  the  assembly  process  of  aligning 
the  transducer  face  to  the  zero  mark  on  the  angle  indicator  dial.   It 
is  estimated  that  horizontal  separations  were  accurate  to  +  0.2  cm. 
The  height  of  the  transducer  above  the  ensonified  sediment  patch  was 
accurate  to  +  0.2  cm  and  the  sediment  surface  within  the  ensonified 
patch  was  level  and  smooth  within  +0.1  cm.   With  a  sound  wavelength  of 
about  1.2  cm,  this  scale  of  roughness  in  bottom  level  is  considered  not 
to  have  a  significant  effect  on  observed  signal  amplitudes.   The  un- 
certainty in  measuring  voltage  amplitude  is  estimated  to  be  2%  of  the 
millivolt  scale  in  Figure  5. 

At  angles  of  incidence  greater  than  70°  ,    small  but  perceptible 
off-axis  energy  from  the  direct  path  was  detected.   Since  the 
difference  between  reflected  and  direct  path  lengths  decreases  with 
greater  angle  of  incidence,  the  arrival  time  difference  between  the  two 
pulses  decreased  so  that  the  leading  edge  of  the  reflected  pulse  over- 
lapped the  tail  of  the  direct  path  off-axis  pulse.   This  resulted  in 
an  estimated  +  5%  error  in  reflected  voltage  amplitudes  above  70°. 

E.   THEORY  OF  MEASUREMENT 

Referring  to  Figure  2,  it  was  assumed  that  the  sound  pressure  ampli- 
tude varied  inversely  as  the  square  of  the  range  from  the  source  and 
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that  the  absorption  in  the  water  was  negligible.   The  assumption  was 
also  made  that  the  sediment  surface  was  sufficiently  smooth  and  homo- 
geneous so  that  the  reflected  rays  behaved  as  if  they  had  originated 
at  the  mirror  image  position  of  the  source. 

The  received  intensity  I,  at  any  range  x  along  the  direct  path,  was 
determined  by  orienting  the  beam  axes  of  the  transducers  toward  each 

other  and  measuring  signal  amplitude  A'.   Intensity  I  is  proportional  to 

2 
(A1)  .   The  reflected  wave  intensity  I  ,  which  is  proportional  to  the 

square  of  the  measured  reflected  wave  amp] itude  A',  is  a  function  of 

the  reflection  coefficient  R,  I  ,  and  the  spherical  spreading  along 

the  total  reflected  path  length  y  =  y  +  y   =  2y  ,  assuming  that  the 

reflection  is  specular.   I   is  the  axial  intensity  at  unit  distance 

from  the  sound  projector,.   Stated  in  equation  form,  this  relation  is 

2        2 
I   =  (A")  -   RI  /y  .   Solving  for  R  and  substituting  for  I   determined 

from  the  direct  path  measurements  gives  the  following: 


2A"y  \  2 
A'x 


Bottom  loss  is  calculated  by 


BL  =  -10  log   R  (dB) 
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IV.   PRESENTATION  OF  DATA 

Tables  I  and  II,  which  list  some  of  the  mass  physical  properties 

of  the  sediment,  are  self-explanatory.   Table  II  records  the  speed  of 

sound  and  density  in  the  sediment  for  three  different  depths  near  the 

surface.   Sound  speed  ratio  c,/c  was  measured  with  a  sound  velocimeter 

1  w 

and  densities  were  measured  using  the  sampling  device  described  by 
Martinek  [5] .   The  data  show  both  density  and  velocity  gradients  in 
the  upper  part  of  the  sediment. 

Table  III  lists  a  series  of  eight  combinations  of  the  five  measured 
sediment  parameters  which  are  used  in  the  theoretical  model.   These 
combinations  are  designed  to  determine  the  model's  sensitivity  to  para- 
meter changes  with  each  combination  given  a  sequence  number  for  refer- 
ence purposes.   Sequence  (1)  is  based  on  Table  I  dynamic  rigidities  at 

3 
a  density  of  1.35  gm/cm  along  with  sound  speed  and  absorption  coeffi- 

3 
cient  at  a  density  of  1.38  gm/cm  .   The  maximum  positive  experimental 

error  from  Table  I  dynamic  rigidities  is  the  change  inserted  in  Sequence 

3 
(2).   The  effect  of  only  a  density  change  from  1.38  gm/cm   to  1.28  gm/ 

3 
cm  is  determined  from  using  Sequence  (3).   Together,  Sequence  (2)  and 

(3)  test  the  model  for  sensitivity  to  density  changes.   Sequence  (4) 

3 

shifts  to  the  Table  I  values  of  dynamic  rigidity  for  1.28  gm/cm  density 

and  in  Sequence  (5)  the  measured  sound  speed  for  a  density  of  1.28  gm/ 

3 
cm  in  Table  III  is  inserted.   Together,  Sequence  (4)  and  (5)  determine 

the  effect  of  only  a  sediment  sound  speed  change.   Also,  Sequence  (5) 
is  based  on  the  measured  data  which  is  assumed  to  be  the  closest  approxi- 
mation to  the  sediment  conditions  during  the  reflection  experiment. 
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Sequence  (6)  and  (7)  are  designed  to  test  the  model  for  large  arbi- 

3 
trary  changes  in  dynamic  rigidity.   The  density  of  1.24  gm/cm   is  an 

assumed  value  of  sediment  density  at  a  point  closer  to  the  surface 

than  that  at  which  a  suitable  sample  could  be  collected.   Sequence  (8) 

tests  the  model  for  a  large  arbitrary  change  in  absorption  coefficient 

relative  to  that  in  Sequence  (5) . 

Tables  IVa  and  IVb  tabulate  the  values  of  R  and  BL  for  the  eight 
sequences  computed  from  Equation  (A-ll)  for  a  range  of  angles  of  inci- 
dence from  14°  to  83°.   From  these  two  tables  sensitivity  of  the 
theoretical  model  is  determined.   The  values  from  Sequences  (1)  and  (5) 
are  plotted  in  Figure  6  for  comparison. 

Table  V  tabulates  the  data  from  the  reflection  experiment.   The 
values  of  BL  in  this  table  are  plotted  in  Figure  6  along  with  calcu- 
lated values  from  Sequences  (1)  and  (5) . 
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TABLE  I 
Physical  Properties  of  Saturated  Water-Kaolin  Sediment 

Cohen  [4]  Sediment  Properties  for  20°C 


Density 

Dynamic  Rigidity 

Pl 
(gm/cm  ) 

-4 
x  10 

2 

(dynes/cm  ) 

-4 
x  10 

2 
(dynes/cm  ) 

1.28 
1.31 
1.35 

26.6  +  3.4 

54.8  +  7.4 

67.9  +14.5 

2.85  +  1.20 

9.68  +  4.36 

16.6   +  6.90 

Martinek  [5]  Sediment  Properties  for  19.9°C 


Frequency 

Density 

Sound  Speed 
** 

Absorption 
Coefficient 

f 
(kHz) 

Pl 

3 
(gm/cm  ) 

°1 
(cm/sec) 

a 
P 

(nepers/cm) 

100 

125 

*    120 

1.38 
1.38 
1.38 

1422  x  102 
1422  x  102 
1422  x  102 

3.79  x  10~3 
4.83  x  10~3 
4.62  x  10~3 

** 


Linear  interpolation  from  100  and  12  5  kHz  data 


c./c   =  0.94 
1  w 
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TABLE  II 
Measured  Sound  Speed  and  Densities 


Average 

Sediment 
Depth 

Average 

Sediment 

Density 

Cl 

—  at  Depth 
c 
w 

** 
Cl 

(cm) 

(gm/cm  ) 

(cm) 

(cm/sec) 

1.55 
(Surface  sample) 

1.28  +  .05 

0.974  at  1.37 

1442  x  102 

3.85 

1.38  +  .05 

0.965  at  3.39 

1428  x  102 

7.85 

1.42  +  .05 

Not  measured 

- 

Measured  at  18.8°C 


**  2 

Corrected  to  20°C,  based  on  c   =  1480  x  10   cm/sec 


w 


TABLE  III 

Combinations  of  Sediment  Parameters 
Used  in  Theoretical  Computations 


Sequence 
No. 

(dynes/ 

2;"x 

cm  ) 

Pl    3 

(gm/cm  ) 

Cl 
(cm/sec) 

a 
P 
(nepers/cm) 

I. 

67.9xl04 

16.6x10 

1.38 

lo422xl05 

-4 
0.462x10 

2. 

82.4xl04 

23.5xl04 

1.38 

1.422xl05 

0.462xl0~4 

30 

82.4x10 

23.5xl04 

1.28 

1.422xl05 

0.462xl0~ 

4. 

26.6xl04 

4 
2.85x10 

1.28 

1.422xl05 

-4 
0.462x10 

5. 

26.6x10 

2.85xl04 

1.28 

1.442xl05 

0.462xl0~4 

6. 

266 

28.5 

1.24 

1.450xl05 

0.462xl0~4 

7. 

8 
26.6x10 

2.85xl08 

1.24 

1.450xl05 

0.462xl0~ 

8. 

26.6x10 

2.85x10 

1.28 

1.442xl05 

0.462xl0~2 
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V.   COMPARISON  OF  DATA 

A  comparison  of  the  bottom  loss  calculations  from  the  reflection 
model  with  the  measured  results  from  the  experiment  reveals  somewhat 
poor  agreement  between  the  two.   The  experimental  limits  on  the  angles 
of  incidence  prevented  measuring  reflection-  coefficients  in  the 
vicinity  of  the  theoretically  predicted  values  of  40°  to  43°.   Figure  6 
shows  a  possible  measured  angle  of  intromission  in  the  vicinity  of  65°, 
but  the  lack  of  experimental  results  below  50°  leaves  open  the  possi- 
bility of  even  higher  measured  bottom  losses  at  the  lower  angles.   An 
additional  comparison  was  made  with  a  liquid-liquid  model  described  in 
Kinsler  and  Frey  [7,  p.  145],  using  only  the  densities  and  sound  speeds 
of  Sequences  (1)  and  (5) ,  which  predicted  an  angle  of  intromission  in 
the  range  of  26°  to  27°.   The  theoretical  magnitudes  of  bottom  loss 
near  the  angle  of  intromission  are  10  to  20  dB  higher  than  the  values 
measured  at  about  65". 

The  eight  combinations  of  parameter  changes  reveal  a  number  of 
facts  regarding  the  theoretical  model  and  its  application  to  the  experi- 
ment performed.   A  summary  of  comparisons  among  sequences  is  contained 
in  Table  VI  on  the  following  page.   The  high  sensitivity  to  density  and 
sound  speed  changes  shown  in  this  table  agrees  with  what  is  expected 
from  the  basic  liquid-liquid  model  [7] .   The  Table  VI  comparisons  also 
reveal  that  dynamic  rigidity  and  absorption  coefficient  affect  bottom 
loss  magnitude  only  near  the  angle  of  intromission.   The  range  of  para- 
meter changes  selected  for  the  sensitivity  tests  had  only  a  small  effect 
on  the  range  of  theoretical  angles  of  intromission,  as  shown  in  Table  VI 
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TABLE  VI 
Comparison  of  Sequences 


Sequence 
Numbers 

1  and  2 


2  and  3 

4  and  5 
6  and  7 

5  and  8 


Parameter 
Change 

Small  change  in 
dynamic  rigidity 


Density 

Sound  Speed 

Large  change  in 
dynamic  rigidity 


Absorption 
Coefficient 


Sensitivity  of 
Equation  (A-ll) 

Sensitive  only  at 
the  angle  of  intro- 
mission 

Very  sensitive 

Very  sensitive 

Sensitive  only  near 
the  angle  of  intro- 
mission 

Low  sensitivity  only 
at  angle  of  intromis- 
sion 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  data  collected  for  one  sediment  condition  at  one  frequency  for 
a  limited  range  of  angles  of  incidence  must  be  considered  as  preliminary 
in  nature.   Three  extensions  of  data  collection  using  the  developed 
experimental  method  are  needed  to  complete  the  experiments.   A  variation 
in  sediment  conditions  is  needed  to  measure  the  effect  on  reflection 
coefficient  from  known  changes  in  density,  sound  speed  and  dynamic 
rigidity.   Apparatus  now  exists  to  measure  all  of  these  parameters.   A 
variation  in  frequency  also  is  needed  to  measure  the  effect  of  changes 
in  absorption  coefficient,  which  has  been  assumed  to  be  the  only  fre- 
quency dependent  sediment  property.   Measurements  over  a  wider  range  of 
angles  of  incidence  are  needed  in  order  to  determine  experimentally  the 
existence  of  an  angle  of  intromission.   In  order  to  increase  the  range 
of  angles  of  incidence  changes  in  the  apparatus  will  be  required.   As 
an  example,  a  deeper  and  larger  tank  will  permit  measurements  at  smaller 
angles  of  incidence. 

The  low  sensitivity  of  R  to  values  of  dynamic  rigidity  could  be  due 
to  the  range  of  dynamic  rigidities  used  for  the  material  not  being  repre- 
sentative.  Because  of  the  fact  that  the  reflection  model  is  quite 
sensitive  to  sound  speed  and  density,  there  exists  the  possibility  that 
a  part  of  the  discrepancy  between  calculated  and  observed  reflection 
coefficients  is  due  to  inadequate  accuracy  for  measured  sediment  physi- 
cal properties. 

The  sediment  densities  and  velocities  measured  during  this  experi- 
ment and  earlier  by  Martinek  [5]  confirm  that  gradients  in  properties 


38 


do  exist  in  the  upper  portion  of  the  sediment.   Although  the  gradients 
may  have  a  significant  influence  on  the  observed  reflection,  this  was 
not  accounted  for  in  the  theoretical  model  used  here  and  time  limita- 
tions did  not  permit  the  adaptation  of  another  model.   One  such  model 
which  does  include  sound  velocity  gradients  in  the  reflecting  medium 
has  been  developed  by  Morris  [8]  based  on  the  theory  of  layered  media 
reflections  described  by  Brekhovskikh  [9] .  .The  gradient  models  must  be 
included  for  experiments  in  which  the  sound  reflection  process  involves 
greater  penetration  of  sound  into  the  sediment  than  was  probably  the 
case  in  the  experiment  reported  here. 
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APPENDIX  A 

Algebraic  steps  to  solve  for  reflection  coefficient  in  Section  II 
are  set  forth  in  this  appendix. 


,.2  ..    .     ,2     2,1/2   -iG      -i( 
(a)   =  x-iy  =  (x  +  y  )     e    =  re 


(A-l) 


Note:   0  not  related  to 


x  E  X^  +  2^ 


y  E  X"^    +    2]l"1 


-    ,    2     2.1/2 
r  =  (x   +  y  ) 


(A- 2) 
(A-3) 
(A- 4) 


a  =  (r)1/2  e"1  9/2  =  (r)1/2(cos  0/2  -  i  sin  6/2) 


The  positive  roots  of  the  following  trigonometric  identities  were 
substituted: 


0  tr>              ,1  +  cos  6/2,1/2    .   .  ,„      ,1  -  cos  6/2,1/2 
cos  6/2  =  +  ( — — '—)         ,    sin  6/2  -  +  ( — ) 


a-  (r)1/2  {[(ijd+i)]^  _    I     x  V2} 

2     r  2     r 


1   r/    xl/2    .  ,    ,1/2, 
[(r+x)     -  i(r-x)    ] 


^" 


Rationalizing  the  inverse  of  the  above  to  obtain  —  yields  the 

a 


following 


,1/2    .,    .1/2 
1    (r+x)     +  i(r-x) 

a  "  rz 

r  V2 
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Substituting  the  real  and  imaginary  parts  into  Equations  (17) 
and  (18)  yields  the  following: 


Re[i)    = 

(r+x)x/" 
r  /2 

lm[i]    = 

(r-x) 
r  \f? 

c.,    = 

r  /2 

/r  +  x 


(A-5) 


a  =  v/2  Trf  ,/r-x 

P       r 


Algebraic  manipulation  of  cn  and  a   then  was  used  to  eliminate 

1      P 


r  and  state  all  relations  in  terms  of  x  and  y. 

2    2r2  2  (7!f  )2(r-x) 


c,   =  ,  a  = 


1    r+x  p       r 


2     2 
r  +  x  -  (-—)  r   =  0 

Cl 


1/  Px2   2   rt 

r-x  -  -(^)   r  =  0 


2x  +  r2[|(i)2__i_]=0 


Cl 


—  -  r—  -  -(i)2!  (A-7) 

2    L   2    4W  J  lA   ' 

Cl 
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Manipulating   further   as   follows: 


(Ci)2(a    )2   =    (^Trf)     (r-x) 
1  p  r+x 


c   a 

,    1  Ps2   _   r-x 

V2TTf    ;        "   r+x 


From 


r      = 


2  2 

x      +   y    , 


x 

1   - 

c   a 
(_LR,  2 

r 

1   + 

c  a 
(_LP)2 

ft2 

r 

=    1   - 

r 

<r 


1- 


K2-nt   '  4nf   ; 


1   +    2     (JLE,2    +    (J^)4 
2'fVf  2Tlf 


ft 

r 


=   + 


c.a 

_J_P 

TTf 


ca 

i  +  (AV 


l27Tf 


and   since   a      and   c      are   greater   than   zero,    the  positive   root   is   chosen, 

IT  "^ 

Combining  x/r  and  y/r  to  obtain  y/x  yields  the  following: 


c,  a 
1  P 

TTf 


c,a 


(A-8) 


1  - 


(_LH)2 


Substituting  the  definitions  of  x,  y  and  r   in  Equations  (A-2) , 
(A-3)  and  (A-4)  shows  that  y '  ,  y"  ,  A1   and  A"   are  related  only  to 
the  groupings  of  known  physical  properties  of  Equations  (A-7)  and  (A-8) 
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Solving  the  interface  equations  by  adding  Equation  (9)  and  (1J ) 
and  substituting  the  result  into  Equation  (10)  gives: 


2a2(V2  "  AV  =  B'l(bl   +  1} 


(b2  -  1)  (b2  -1)B' 
P2  fc)   (A-a  +  A-2)  -  Vl[   •  +  2b1B'1 


Substituting 


Z  = 


B' 


A'   +  A". 


A'   -  A" 

A  2   A  2 


2a2  (A'2  -  AV 
(b2x  +  1) 


2Vik 


2  „  2 
P2W  (b    +  1) 


2       2 
2a, 


+  2b. 


(A-9) 


Letting  Z  equal  the  right  hand  side  of  Equation  (A-9)  and  defining 
intensity  reflection  coefficient 

A"2 

-r-~-  (A-10) 


such  that 


A'2+A"2    1+A»2/A'2    1+/J 


A'2  -  A"2    1  -  A"2/A'2    1  -  /r 


=  Z 


then 


R  = 


Z-l 
Z+l 


(A-ll) 
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All  quantities  on  the  right  hand  side  of  Equation  (A-ll)  are  physi- 
cal properties  of  the  sediment  or  water  and  Z  is  a  complex  number. 
Formulating  a   and  b  as  complex  numbers  gives  the  following: 


/    w   .2 


H  and  J   are   defined   as    follows: 


Cw      .2 


pi  to"e7>     a'i  +  2yV 


H      =        ±— -    1  (A-13) 

(A^  +  2y'1)    +(X"1  +  2y"1) 


w       ,  2 


Pi    <ii^TT>       (A"l   +   ^Y 


j   =  ± (A- 14) 

(X'j  +  2y'1)    +(A"1  +  2y"1) 


.         w         2 

pi  (ii^eT5 

b    ,    = 1  =   F  +    iG  (A-15) 


1 


v\  -  iy-L 


F  and   G   are  defined   as   follows: 


c  o 

I        ,         w       ^ 2 

piy  i  (iirTe-) 

F   =        5 ^T      "    1  (A_16) 

(yy     +  (y'^) 


„        .         w .  2 

Kl      1      sin   o. 

l 


G 


(y^)2  +  (y^)2 
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At  this  point  the  two  Lame  constants  A'  and  A"  are  contained  only 
in  the  terms  H  and  J.  Equations  (A-13)  and  (A-14)  are  then  restated  in 
terms  of  x  and  y  as  defined  earlier  in  Equations  (A-2)  and  (A-3). 


pix  'iLTeT'              p.x  c2 
■"  2 T~    -  l=— 1 

x   +  y  r 


J  = 


P1Y  (sTn\:)2     pv  c2 
1        1 


2     2 
x   +  y 


2 


j  =  p  c2(-_)  A  -  A  (h  +  l) 

1      2   x      x 


Substituting  from  Equations  (A-7)  and  (A-8)  then  yields  the  final 
expressions  for  H  and  J. 


J  = 


/   w  .2 
P.  (~ — s~ ) 

1  sin  6. 


c  a 
TTf 


i(J-)2 


V27Tf 


-  1 


H  +  1 


(A-18) 


(A-19) 
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